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Abstract. The discrepancy between the results of covariant and non- 
covariant one-loop calculations for higher-spin fields in quantum cosmol- 
ogy is analyzed. A detailed mode-by-mode study of perturbative quantum 
gravity about a flat Euclidean background bounded by two concentric 3- 
spheres, including non-physical degrees of freedom and ghost modes, leads 
to one-loop amplitudes in agreement with the covariant Schwinger-DeWitt 
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method. This calculation provides the generalization of a previous analysis 
of fermionic fields and electromagnetic fields at one-loop about flat Euclidean 
backgrounds admitting a well-defined 3+1 decomposition. 
PACS numbers: 03.70. +k, 04.60. +n, 98.80.Hw 

I. INTRODUCTION 

In recent papers devoted to one-loop calculations in quantum gravity and 
quantum cosmology [1-11], discrepancies were found between the covariant 
Schwinger-DeWitt method [12], where the scaling factor of one-loop ampli- 
tudes coincides with the A 2 coefficient in the heat-kernel expansion, and 
zeta-function regularization [13-14] with the corresponding mode-by-mode 
analysis of quantized fields. 

In Ref. [1] it was shown that the C(0) value for gravitons calculated on 
the Riemannian de Sitter 4-sphere in terms of physical degrees of freedom 
[15] disagrees with that obtained in covariant formalism [16]. In Refs. [2-11] 
it was shown that the results of analogous calculations for gauge fields on 
manifolds with boundaries, as well as earlier results of Refs. [17-18], are in 
disagreement with the results of covariant calculations using the modified 
Schwinger-DeWitt formulae for manifolds with boundaries [19]. 

In Ref. [20] the attempt of investigating the reason of the discrepancy 
described above was made. Since for spin-| fields there is no gauge freedom, 
different results cannot be due to inequivalent quantization techniques, but 
may have an entirely geometrical origin. The hypothesis was then put for- 
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ward that the reason of this discrepancy consists in the inappropriate use of 
3+1 decomposition on the manifolds where such a split is ill-defined (appar- 
ently the 3+1 decomposition is ill-defined on the full 4-dimensional sphere 
or on the part of 4-sphere or 4-dimensional flat space bounded by a single 3- 
sphere, since the vector field matching the normal to the boundary is singular 
at the origin). The £(0) value for massless spin-1/2 fields on the manifold 
representing the part of flat Euclidean 4-dimensional space bounded by two 
concentric 3-spheres (such a manifold admits a well-defined 3+1 split) was 
calculated. It was shown that the result coincides with the covariant one 
obtained within the framework of the Schwinger-DeWitt method [3,19]. 

On evaluating £(0) for vacuum Maxwell theory subject to the Coulomb 
gauge before quantization, the one-loop result is in disagreement with the 
covariant one, even in the two-boundary case, where one has a well-defined 
3+1 split [21]. The possible reason is that, on non-trivial backgrounds with 
boundaries, contributions of ghosts and non-physical degrees of freedom do 
not cancel each other. In Ref. [21], continuing the study of the quantization of 
the electromagnetic field on manifolds with boundaries appearing in Ref. [22], 
this suggestion was checked. It was shown that, if one deals with non-physical 
degrees of freedom and ghost modes within the Faddeev-Popov formalism, 
and if one studies a flat Euclidean background bounded by two concentric 
3-spheres, the one-loop evaluation of quantum amplitudes agrees with the 
covariant result, confirming in such a way our hypotheses. A technique for 
the disentanglement of coupled gauge modes was also described and applied 



Gravitons in one-loop quantum cosmology 



4 



for the first time. 

Here, we calculate £(0) for gravitons on the part of flat 4- dimensional Eu- 
clidean space bounded by two concentric 3-spheres, taking into account the 
contributions of non-physical modes and ghosts. We show that in this case, 
which is more complicated from the technical point of view than the electro- 
magnetic analysis of Refs. [18,21-22], one finds again agreement between the 
results of covariant and non-covariant calculations. 

In our calculations we use the version of the generalized ^-function tech- 
nique elaborated in [8-10]. The main ideas are as follows. Zeta-functions 
are traces of complex powers of elliptic, self-adjoint, positive-definite differ- 
ential operators. They have an analytic continuation to the complex plane as 
meromorphic functions with simple poles. Remarkably, they are regular at 
the origin together with their first derivative, and this enables one to define 
and compute the determinants of the corresponding operators and one-loop 
quantum amplitudes. With our notation, one writes f n (M 2 ) for the function 
occurring in the equation obeyed by the eigenvalues by virtue of boundary 
conditions, and d(n) for the degeneracy of the eigenvalues. One then defines 
the function 

oo 

I(M 2 ,s) = £ d(n)n- 2s log f n (M 2 ). (1.1) 

n=no 

Such a function has a unique analytic continuation to the whole complex-s 
plane as a meromorphic function, i.e. 

I(M 2 , s) = /polc(M2) + I R (M 2 ) + O(s). (1.2) 
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The £(0) value is then obtained as [8-10] 



C(0) = /log + ipolc(oo) - ipole(O), 



(1.3) 



where ii og = l£ is the coefficient of logM from I(M 2 , s) as M — > oo, and 
7 po ie(M 2 ) is the residue at s = 0. Remarkably, 7i og and 7 po i e (oo) are obtained 
from the uniform asymptotic expansions of modified Bessel functions as their 
order tends to oo and M — > oo, whereas / po i e (0) is obtained from the limiting 
behaviour of such Bessel functions as M — > 0. 

In Sec. II we write down the equations for basis functions for physi- 
cal and non-physical degrees of freedom of the gravitational field with the 
corresponding ghost modes, and we find their solutions on the chosen back- 
ground. In Sec. Ill we consider the boundary conditions and compute all 
contributions to the full £(0) value. Results and concluding remarks are pre- 
sented in Sec. IV. The forms of the differential operators acting on metric 
perturbations are given in the Appendix. 



II. EQUATIONS FOR BASIS FUNCTIONS 
AND THEIR SOLUTIONS 



For the reasons described in the introduction, we study pure gravity at 
one-loop about a flat Euclidean background with two concentric 3-sphere 
boundaries. Our approach to quantization follows the Faddeev-Popov for- 
malism (cf. [22]). Hence we deal with quantum amplitudes of the form 
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With our notation, C is the set of all Riemannian 4-geometries matching 
the boundary data, [i\ is a suitable measure on the space of metrics, /i 2 is 
a suitable measure for ghosts, § v is an arbitrary gauge- averaging functional, 
and the total Euclidean action reads (in c = 1 units) 

r * = '* + lib I (4,fl v^ ** + L Tr K ^~ q * x 

+ ^Ll^^^ < 21 > 

Of course, K is the extrinsic-curvature tensor of the boundary, q is the in- 
duced 3-metric of dM, and a is a positive dimensionless parameter. The 
ghost action I g h depends on the specific form of <Ev Denoting by h^ v the 
perturbation around the background 4- metric g^ u , one thus finds equations 
of motion of the kind 

where □* is the 4-dimensional elliptic operator corresponding to the form 
of one is working with. Here we choose the de Donder gauge-averaging 
functional 

= v ( V - \gj) » 

where V is covariant differentiation with respect to g^ u , and h = g*" V- 
The corresponding U^ DD operator is the one obtained by analytic continu- 
ation of the standard D'Alembert operator, hereafter denoted by □. The 
resulting eigenvalue equation is (see Appendix) 

Uh uv + \h uv = 0. 
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Now we can make the 3+1 decomposition of our background 4-geometry 
and expand h 0Q , h 0i and in hyperspherical harmonics as 



n=l 



h 00 (x,T) = Y / a n (r)Q^(x), 

oo 

h i(x,T) = 



n=2 



b n (r) ViQ 2 {n) ^ +c n (r)St ) (x) 



(n 2 - 1) 



(2.2) 
(2.3) 



hij(x : r) = J2 d n( T ) y (n 2 -l) 



3 



+ E fn(r)(v4 n \x) + V,-5f >(*)) + k n (r)G\f(x) 



n=3 



n=l 6 

(2.4) 



Here Q^ n \x), Sl n \x) and G\? (x) are scalar, transverse vector and transverse- 
traceless tensor hyperspherical harmonics respectively, on a unit 3-sphere 
with metric q,-. Their properties are described in Refs. [17-18,23]. 

The insertion of the expansions (2.2)-(2.4) into Eq. (2.1) leads to the 
following system of equations (decoupled modes will be treated separately): 



A n a n (r) + B n b n (r) + C n e n {r) = 0, 

D n b n {r) + E n a n (r) + F n d n (r) + G n e n (r) = 0, 
L n d n {r) + M n b n {r) = 0, 
N n e n (r) + P n b n (r) + Q n a n (r) = 0, 
H n c n {r) + K n fn{r) = 0, 
R n fn(r) + S n c n (r) = 0, 
f n k n (r) = 0. 



(2.5 

(2.6 
(2.7 
(2.8 
(2.9 
(2.10 
(2.11 
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Since our background is flat, after setting a = 1 in (2.1) the operators ap- 
pearing in Eqs. (2.5)-(2.11) take the form (for all integer n > 3) 

2 . d 2 3d (n 2 + 5) 



D n = — + -—- ± 5-^ + A 



d 2 

dr 2 + 


3 d 


(n 2 + 5) 


t dr 


r 2 




B n = 


4 
: r 3 "' 




= 


2 

r 1 ' 


a 

dr 2 + 


1 a 


(n z + 4) 


t drr 


r 2 


E n 


2 

s 7< 


^ 2 -l), 


F n 


~~ 3 


r? 2 - 4) 
r 3 ' 


G n = 


o 
z 

3 


(n 2 1 "I 
T 6 


d 2 

dr 2 + 


1 d 


in 2 + h) 


r dr 


r 2 


K 


2 


'n 2 - 4) 


d 2 


1 d 


(n 2 - 5) 


dr 2 


r dr 


r 2 




M„ = 


_ 4 
r 


d 2 


1 d 


(n 2 + 1) 


dr 2 


r dr 


r 2 






4 
r 




= 


= 6, 


d 2 


1 d 


(n 2 - 4) 


dr 2 


r dr 


r 2 



# n = — + -— - ^ —l + A 



L n = -rr; : n 1" A 



Rn = -rr; ; n 1" A 
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S n = —■ 

T 



~ d 2 1 d (n 2 - 1) , , 

Inserting the operator T n from Eq. (2.12) into Eq. (2.11) we can easily find 
the basis function describing the transverse-traceless symmetric tensor har- 
monics which usually are treated as physical degrees of freedom [1,15,17,23] 

k n (r) = «ir/ n (Mr) +a 2 rK n (Mr), n = 3, . . . (2.13) 

where M = \/—X and / and K are modified Bessel functions. 

However, the equations (2.5)-(2.8) for scalar-type gravitational perturba- 
tions lead to a rather complicated entangled system as well as Eqs. (2.9)- 
(2.10), describing vector perturbations. In Ref. [21], where we have studied 
the analogous problem for the electromagnetic field, a method was used to 
decouple a similar entangled system for normal and longitudinal components 
of the 4-vector potential. The idea is that one can diagonalize a 2 x 2 op- 
erator matrix after multiplying it by two functional matrices. In some cases 
one can choose these functional matrices in such a way that the transformed 
operator matrix is diagonal and the corresponding differential equations for 
basis functions are decoupled. However, in the case of scalar-type gravita- 
tional perturbations we have a 4 x 4 operator matrix. To diagonalize such 
a matrix it is necessary to solve a system of 24 second-order algebraic equa- 
tions with 24 variables. This problem seems a rather cumbersome one and 
we thus use another method. For this purpose, we assume that the solution 
of the system of equations (2.5)-(2.8) is some set of modified Bessel functions 
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with unknown index v. Let us look for a solution of this system in the form 

aM^^^l, (2.14) 



b n (r) = foW v (MT), (2.15) 

d n (r)=(3 3 TW u (Mr), (2.16) 

e n {r)=(3 A rW v {Mr). (2.17) 

Here, W v is a linear combination of modified Bessel functions I v and K v 
obeying the Bessel equation 

Now, inserting the functions (2.14)-(2.17) and the corresponding operators 
from Eq. (2.12) into the system of equations (2.5)-(2.8), and taking into ac- 
count the Bessel equation (2.18), one finds the following system of equations 

for /?i , /? 2 , /?3 and /3 4 : 

{v 2 -n 2 - 6)/3i + 4/5 2 + 2/3 4 = 0, 

6(n 2 - 1)A + 3(z/ 2 - n 2 - 4)/? 2 + 4(n 2 - 4)/3 3 - 2(n 2 - l)/3 4 = 0, 
4/3 2 + (z/ 2 -n 2 + 4)/? 3 = 0, 

6/?i - 4/? 2 + (z/ 2 - n 2 - 2)04 = 0. (2.19) 

The condition for the existence of nontrivial solutions of the system (2.19) is 
the vanishing of its determinant, i.e. 

(z/ 2 - n 2 ) 2 \{u 2 - n 2 ) 2 - 8(z/ 2 - n 2 ) - 16(n 2 - 1)1 = 0. (2.20) 
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The roots of Eq. (2.20) are 

v 2 = n 2 , v 2 = (n - 2) 2 , v 2 = (n + 2) 2 . 

The positive values of v provide the orders of modified Bessel functions. Now 
we can write down the /3's corresponding to different values for z/'s. For v = n 
one has 

ft = 3A, /3 2 = /?s = 0, (2.21) 

or 

A = 0, (h = -fo,P4 = -2fo. (2.22) 

For v = n — 2 one has 

02 = (n + A = r^4k, 0* = -01- (2-23) 
Last, for z/ = n + 2 one has 

(n — 1) 

2 = -{n - 1)0,, 3 = T^j0i, 04 = -01- (2.24) 

Having the Eqs. (2.21)-(2.24) we can get the basis functions for scalar-type 
gravitational perturbations (2.14)-(2.17) 

a n (r) = - (7i4(Mr) + 7 3 / n _ 2 (Mr) + 74 / n+2 (MT) 

+5 1 K n (Mr) + 5 3 X n - 2 (Mr) + 5 4 K n+2 (Mr)), (2.25) 



6n(r) = 7 2 4(Mr) + (n + l)7 3 / n _ 2 (Mr) 
-(n - l) 7 4/n +2 (Mr) + 5 2 K n (Mr) 

+ (n + l)5 3 K n _ 2 (Mr) - (n - l)5 4 K n+2 (Mr), (2.26) 
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d n (r) = r (- l2 I n {Mr) + |^±il 73 /„_ 2 (Mr) 

+|^74/n + 2(Mr) - 5 2 K n (Mr) 
(n + 2) 

' " L) 5 3 K n _ 2 (Mr) + ^—^8 4 K n+2 {Mr)\ (2.27) 



(ra-2) ° ™ 7 (ra + 2) 

e n (r) = r(3 7l / n (Mr) - 2 72 / n (Mr) - 7 3 /„_ 2 (Mr) 
-7 4 4 +2 (Mr) + 3<fitf n (M<r) - 25 2 K n (Mr) 
-5 3 K n . 2 (Mr) - 5 4 K n+2 (Mr)). (2.28) 

We can find the basis functions for vector-like gravitational perturbations 
in a similar way. Let us suppose that 

c n (r) = e 1 W u (Mr) (2.29) 

and 

f n (r)=e 2 rW u (Mr). (2.30) 
Inserting (2.29)-(2.30) into Eqs. (2.9)-(2.10) one has the system 

{y 2 - n 2 - 5)ei + 2(n 2 - 4)e 2 = 0, 

2e t + (u 2 -n 2 + 3)e 2 = 0. (2.31) 
The determinant of the system (2.31) is 

{y 2 - n 2 ) 2 - 2{y 2 - n 2 ) - An 2 + 1 

and its positive roots are nil. For v — n + 1 one has 

1 

" 2 " ~(^f 1 
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and for v = n — 1 one has 

1 

(ra-2) 

and correspondingly the functions (2.29)-(2.30) take the form 

(Mr) + £ 2 / n _i(Mr) + Vl K n+1 (Mr) + 772 (Mr), (2.32) 

/n(r) ^ + e 2 In-i{Mr) 

\ [n + 2) (n — 2) 

-0TT2)" lA '"« (Mr) + (^2)" 2A '"- (Mt) ) ' (2 ' 33) 

We have also to find the basis functions for ghosts. The eigenvalue equa- 
tions for ghosts in the de Donder gauge have the form (see Appendix) 

□y^ + Ay^ = 

and the corresponding fields can be expanded on a family of 3-spheres as 

00 

Mx,t) = 'EUt)QW(x), (2.34) 



n=l 



(Pi(x,r) = 



n=2 



r , ViQ (n) (x) , , w (n), x 

m 4 r ) ( w 2_x) +P"( r )^ ( x ) 



(2.35) 



The functions /„(r),m n (r) and p n (r) can be found similarly to those for 
harmonics of gravitational perturbations. They have the form 

l n (r) = ^(K 1 I n+1 (Mr) + K 2 I n . 1 (Mr)+6 1 K n+1 (Mr)+e 2 K n . 1 (Mr)), (2.36) 

m n (r) = -(n - l^/^Mr) + (n + 1)k 2 4-i(Mt) 

-(n - l)^X n+1 (Mr) + {n + lj^^^Mr), (2.37) 

p n (r) = m n {Mr) + pK n {Mr). (2.38) 
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III. BOUNDARY CONDITIONS 



We use the boundary conditions for linearized gravity studied in Ref. 
[2]. Assuming that the spatial components hij of perturbations of the gravi- 
tational field, and the normal component ipo of the ghost field vanish on the 
boundary, the authors of Ref. [2] obtained boundary conditions for gauge- 
invariant amplitudes in the form 

hij\dM — KoldM — Vo\dM — 0; (3-1) 
(^ + ^oo-^^)|^ = 0, (3.2) 

(^T " \dM = 0. (3.3) 

Inserting the expansions (2.2)-(2.4) and (2.34)-(2.35) into (3.1)-(3.3), one 
finds the following boundary conditions on the basis functions: 

( da n {r) 6a w (r) I dejj)\ 

{-^- + ~r~ ~ v-^r ) ldM = °' (3 - 4) 

b n (r)\ dM = 0, (3.5) 

c n (r)\dM = 0, (3.6) 

d n (i~)\dM = Q, (3-7) 

e n (r)\ dM = 0, (3.8) 

f n (T)\ dM = 0, (3-9) 
k n (T)\ dM = 0, (3.10) 
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(3.11) 

0, (3.12) 

0. (3.13) 

Now we can calculate £(0) for the gravitational field on the background 
which represents the part of 4-dimensional flat Euclidean space bounded by 
two concentric 3-spheres of radii r + and r_ respectively. Let us begin with 
the contribution of the transverse-traceless harmonics of the gravitational 
field to C(0). On inserting (2.13) into the boundary condition (3.10) one gets 
the system of equations 

aj~ + a 2 K~ = 0, 

ai I+ + a 2 K+ = 0, (3.14) 

where 

I- = / n (Mr_),7+ ee I n (Mr + ),K- ee K n {MrJ),K+ ee K n (Mr + ). 

The condition for the existence of non-trivial solutions of the system (3.14) 
is the vanishing of its determinant, i.e. 

I-K+-I+K~ = 0. (3.15) 

Eq. (3.15) provides the eigenvalue condition which we need to apply the 
version of ^-function technique described in [8-10]. It is obvious that we can 
neglect the first term on the left-hand side of Eq. (3.15). Then using the 
power-series expansion for I and K functions [9] one can show that 7 po ie(0) 



ln(r)\dM = 0, 

dm n (r) 2 „ s \ 
dr 

( dp n (r) 2 



-m n {T 9M = 



dr ~T Pn[T) V 9M 
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is equal to zero, while using the uniform asymptotic expansions of modified 
Bessel functions one can show that 7 po i e (oo) vanishes as well (cf. [21]). Thus, 
we only have to calculate I\ og . Looking at the uniform asymptotic expansions 
for modified Bessel functions one can easily see that the coefficient of log M in 
the determinant (3.15) is equal to (—1). Bearing in mind that the degeneracy 
for tensor transverse-traceless harmonics is 2(n 2 — 4) [17,23], one gets 

oo 

^iog = -4)(-l) = -C*(-2) + 4C«(0) + (-3) = -5, 

rz=3 

where Cr( s ) is the usual Riemann ^-function, and we have used its well-known 
values Cr(— 2) = and (r(0) = —1/2. Thus, now we can write 

C(0)tensor = "5. (3.16) 

Let us now evaluate the contribution from vector harmonics. Inserting 
Eqs. (2.32)-(2.33) into the boundary conditions (3.6) and (3.9) one has the 
system of equations 

+ s 2 lti + + mK-i = °> 

~ £l (^+2) + 62 (^2) " Vl (^T2) + V2 (^2) " °' 

7"+ T + K + K+ 
— — + r — Tji — ^ + ^27 = 0. 3.17 

(n + 2) in -2) ' (n + 2) (n - 2) V ; 

The vanishing of the determinant of the system (3.17) is the eigenvalue 
condition for vector harmonics. Taking into account only the dominant terms 
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of this determinant (i.e. the terms including I + and K ), one can write the 
eigenvalue condition in the form 

7 n Vn-l^n + l^n-l = 0, (3.18) 

where we have omitted the unessential common multiplier which does not 
depend on M. One finds that just as in the case of tensor perturbations 

ipole(oo) = Jpolc(O) = 

and we have to calculate 7i og . It is easy to see from uniform asymptotic 
expansions of modified Bessel functions [24] that 

oo 

^iog = E(™ 2 - !)(-2) = -2Cfl("2) + 2Cr(0) + 6 = 5. 

71=3 

Here, 2(n 2 — 1) is the degeneracy of vector harmonics [18,23]. Thus, we have 

C(0)vcctor = 5. (3.19) 

Further to the entangled modes c n (r) and f n (j) where n — 3, ... we have 

also the decoupled mode c 2 (r) obeying the equation 

^ c 2(r) 1 dc 2 {r) 9c 2 (r) 2 , , 
~ll^ + r—T —-Mc 2 (t)=0, 

whose solution can be written as 

C 2 (T)=8h(MT)+< n Kz(MT). 

This function satisfies Dirichlet boundary conditions on the 3-sphere bound- 
aries and the corresponding eigenvalue condition is 

j 3 + k 3 - = 0. 
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Its /jog and correspondingly £(0) is 



vector decoupled 



-3. 



(3.20) 



The most complicated task is the calculation of the contribution of scalar 
harmonics, since in this case we deal with the system of four entangled equa- 
tions for four basis functions a n {r) , b n (r) , d n {r) and e n (r). Inserting the 
expressions for these functions from Eqs. (2.25)-(2.28) into the boundary 
conditions (3.4), (3.5), (3.7) and (3.8) one finds the following system of equa- 
tions at each boundary: 



271(4* - Mt ± I'±) + 2 72 (/± + Mt ± I'±) + 73(6/^2 + ™t ± I'±_ 2 ) 
+74(6/^2 + 2Mr ± + 25 1 (K± - Mt ± K>±) + 26 2 (K± + Mt ± K>±) 
+5 3 (6^_2 + 2Mt ± K>±_ 2 ) + h(6K± +2 + 2Mr ± K' n ± +2 ) = 0, 
72 /* + in + 1)73^-2 - (n - l)74^ +2 
+5 2 K± + (n+ 1)5 3 K±_2 -in- l)5 4 K± +2 = 0, 



The determinant of this system of equations is rather a cumbersome one, 
and we only write down its dominant part, which reads 




3 7 l/^ - 2 72 /± - 73 J±_ 2 - 74^+2 

+35 1 K± - 25 2 K± - 6 3 K±_ 2 - 5 4 K± +2 = 0. 



(3.21) 
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- 6(n ^Mr + (C 2 ) + 4 + + 2 (4 + ) 2 

mt + (/; +2 ) + iu (4 + ) 2 l 



(n + 2) 
6(n 2 - r 
" (n-2) 



6(n 2 - 1) 



2 



n + 2) 

6l "'" L) ^_« +2 )-X-_ 2 (X-) 2 ] (3.22) 



(n-2) 

Inserting into Eq. (3.22) the expressions for uniform asymptotic expansions 
of modified Bessel functions [24] we see that in the limit M — > oo our deter- 
minant becomes 

144n 2 (n 2 - l) 2 
(n 2 - 4) 2 

and after taking the logarithm and expanding it in inverse powers of n we 
only get vanishing contributions to ipoie, and hence 

^poie(oo) = 0. 

Using the power series for modified Bessel functions, in an analogous way we 
see that in the M — > limit the determinant (3.22) becomes 

_144(n 2 -l) 3 (n 2 -16) 
n 2 (n 2 -4) 3 

This leads to vanishing contributions to 7 po i e and hence 

/pole(0) = 0. 

The calculation of 7i og is straightforward and yields 
/io g = Ey(-2) = -Ci ? (-2) + 5 = 5. 

n=3 
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Thus, we can write 

C(0)scalar = 5. (3.23) 

Now we evaluate the contribution of partially decoupled modes. When 
n = 2 we have the system of equations for a 2 (r), 5 2 (r) and e 2 (r), while the 
mode g? 2 (t) does not exist. Even without writing down explicitly the corre- 
sponding determinant, one can show that its dominant part is proportional 
to the product of two functions /, two functions K and one I' multiplied by 
M and one K' multiplied by M too. Hence their contribution to 7i og , which 
in the case of finite number of degrees of freedom coincides with £(0) (see 
Ref. [8]), can be easily obtained knowing the uniform asymptotic expansions 
of I and K and extracting from them the terms proportional to logM in 
logarithms of these functions. Thus 

/i og = n 2 /2 x (-1)U =2 = -2 

and correspondingly 

C(0)scalarn=2 = -2. (3.24) 

When n — 1 only the modes gh(t) and ei(r) survive, and in an analogous 
way one can show that 

C(0)scalar„=l = 0. (3.25) 

The only thing which we have to calculate now is the contribution of 
ghost fields to the full C(0). Inserting the functions / n (r),m„(r) and p n {r) 
from Eqs. (2.36)-(2.38) into the boundary conditions (3.11)-(3.13) we can 
get the determinants for the corresponding system of equations. From these 
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determinants we can obtain, by the method described above, the following 
contributions to C(0): 

C (0) ghost vector = 2' (3.26) 
C (0) ghost scalar = 0, (3.27) 
C(0) ghost decoupled = — ^" (3.28) 

Of course, the ghost contributions (3.26)-(3.28) to the full £(0) should be 
multiplied by —2 [2,21], but this does not affect our result. Thus, by virtue 
of (3.16), (3.19), (3.20), (3.23)-(3.28) the full C(0) value is given by 

C(0)totai = 0. (3.29) 

This result obtained by using the ^-function technique coincides with that 
obtained by using the covariant Schwinger-DeWitt technique for manifolds 
with boundaries [2,19]. In fact, on using the covariant technique on the part 
of flat Euclidean space bounded by two concentric 3-spheres we see that the 
volume contribution to the A 2 Schwinger-DeWitt coefficient vanishes, while 
the surface contributions from the two boundaries cancel each other. 

IV. CONCLUSIONS 

We have shown that, in quantum cosmology, one can obtain unambigu- 
ous calculations of one-loop amplitudes providing one studies flat Euclidean 
backgrounds bounded by two concentric 3-spheres, and providing one takes 
into account non-physical degrees of freedom and ghost modes. One then 
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finds that the covariant Schwinger-DeWitt technique, and zeta-function reg- 
ularization relying on a mode-by-mode analysis of quantized fields, are in 
agreement. 

Moreover, when the background 4-geometry has boundaries, one finds it 
is no longer true that ghost modes cancel the contribution of non-physical 
degrees of freedom [21]. To preserve gauge invariance one has thus to deal 
with physical, non-physical and ghost modes. All these properties have been 
here proved in the case of pure gravity. 

In Ref. [20] it was shown that discrepancies do not occur when massless 
spin-| fields are studied at one-loop about Riemannian 4-geometries with 
two 3-sphere boundaries. In Ref. [21] we have shown that, on studying 
vacuum Maxwell theory within the Faddeev-Popov formalism, discrepancies 
in one-loop calculations are again eliminated providing one takes a flat Eu- 
clidean background bounded by two concentric 3-spheres. Here, we have 
proved that the same property holds in the more involved case of linearized 
gravity. Interestingly, in our paper the contribution of physical degrees of 
freedom, i.e. transverse-traceless gravitons, is corrected by the contribution 
of non-physical degrees of freedom, while the contribution of ghosts on the 
background under consideration vanishes. 

It therefore seems that, on considering Riemannian 4-geometries with two 
boundaries, covariant or mode-by-mode descriptions of quantum amplitudes 
are both legitimate, providing one takes into account non-physical degrees of 
freedom and ghost fields. By contrast, on studying Lorentzian 4-geometries, 
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the 3+1 decomposition and the extraction of physical degrees of freedom is 
still valid, providing one can make sense of the corresponding path integral. 
The Lorentzian regime is in turn more relevant for the description of the 
time-evolution of the physical universe studied within the framework of the 
Hartle-Hawking program [25-28]. 

Last, but not least, we would like to emphasize that the mode-by-mode 
quantization program may shed new light on modern quantum field theory. 
We are currently investigating the relation between the I(M 2 ,s) function 
defined in Eq. (1.1) [8-11], and the zeta-function. This would enable one to 
relate gauge invariance of quantum amplitudes in the presence of boundaries 
to the invariance under homotopy of the residue of a meromorphic function. 
In physical language, this happens since a change of the gauge-averaging 
functional leads to a smooth variation of the corresponding matrix of elliptic 
self-adjoint operators. The residues at the origin of the meromorphic func- 
tions occurring in this analysis may be studied by a suitable generalization of 
the Atiyah-Patodi-Singer theory of Riemannian 4-geometries with boundary 
[29]. Although it is unclear whether such a research program can be com- 
pleted, and applied to linearized gravity, it seems to provide a very exciting 
and deep vision of gauge invariance in quantum field theory. 

* * * 

After submitting our paper, we became aware of Ref. [30], where the 
asymptotic expansion of the heat kernel for elliptic operators on Riemannian 
4-manifolds with boundary is studied to improve the analysis in Ref. [19]. 
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In Ref. [31] Moss and Poletti, using the new results in Ref. [30], have re- 
calculated the conformal anomalies on Einstein spaces with boundary and, 
in the 1-boundary case (i.e. the disk), they have found agreement with 
the results obtained in Refs. [5,11,20] for spin-| fields, and in Ref. [21] 
for vacuum Maxwell theory. However, it should be emphasized that the 
agreement between the covariant Schwinger-DeWitt and the non-covariant 
mode-by-mode calculations of the covariant Faddeev-Popov path integral is 
achieved on the disk in the Lorentz gauge only. In this gauge, the second- 
order differential operator for the electromagnetic field is covariant and does 
not depend on the choice of 3+1 decomposition. For a more general class 
of relativistic gauges [21-22], the mode-by-mode analysis of 1-loop quantum 
amplitudes is gauge-dependent on the disk, while in the 2-boundary case (i.e. 
the ring), the £(0) value is gauge-independent. 

As far as we can see, this last remaining discrepancy seems to point out to 
serious limitations of the quantum theory when the background 4-geometry 
does not admit a well-defined 3+1 decomposition. Nevertheless, in the light 
of the recent literature [30-31], no conclusive argument exists which proves 
that the gauge-invariance problem can only be addressed in the 2-boundary 
case. It also appears interesting to compare our mixed boundary conditions 
for linearized gravity with the boundary conditions studied in Ref. [32]. 
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APPENDIX 



In Sec. II, the eigenvalue equations (2.5)-(2.11) are obtained out of the 
□ operator, which is the elliptic operator defined by 

n = g ^V,V v . (A.l) 

As stated in Sec. II, covariant differentiation is performed with respect 
to the flat 4-metric g, in the spherical local coordinates suitable for the 
description of flat Euclidean 4-space bounded by two concentric 3-spheres. 
The corresponding perturbation of g is denoted by h, and r is the Euclidean- 
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time coordinate. Hence one finds 

d 2 h 00 3 dh 00 1 \i 4 {i 6 2 .. 

D/i oo = -^2- + --^r + ^oo|i - ^ V - ^oo + (A.2) 

„, d 2 h 0k ldh ok 1 |j 7 2 |j 2 , , 

ofe = + + ^Kk\i - ~2 h M - ^ h ik + - V*> ( A - 3 ) 

d 2 /i i7 - 19k- 1 , | fc 2 , 2 / , \ 2 , IK AX 

nh ^ = -^T---Q^ + -2^\k -- 2 h l3 + -[h Mj +h m ) + - 9lJ h^ (A.4) 

where a vertical stroke denotes 3-dimensional covariant differentiation on a 

3-sphere of unit radius [17-18]. In a similar way, the ghost operators acting 
on (2.34)-(2.35) are found to be (cf. [21]) 

d 2 (f 3 d(p 1 u 3 2 i fc /*x 

= -^r + + ^ - ~2 ^ ~ ~z W , (A.5) 

#Vi 1 dcpi 1 | fe 2 2 

D ^ = ^ + -^ + ^^|fc -~2^ + -^- (A.6) 
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